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1. INTRODUCTION 


1.1. STATEMENT OF THE PROBLEM 


- -The flow of fluid about, a solid. body can be separated 
into two regimes: a thin layer of fluid adjacent to the body 
surface where friction is the dominant phenomenon , and the 
region outside the thin layer where the fluid can be treated 
as inviscid. The thin layer adjacent to the body surface is 
called the velocity boundary layer. The thickness of the 
boundary layer ■,& y increases along the body in a downstream 
direction.- Inside the boundary layer the flow is retarded 
because of friction with the solid surface and it forms a 
velocity profile that varies smoothly from zero at the wall 
to- the undisturbed free stream velocity U Q . The reduction of 
skin friction drag between a solid surface and a moving fluid 
can greatly improve the performance of many vehicles or fluid 
machinery. 

Boundary layer development over a flat plate can be 
separated into three regions. A laminar region starts from 
the flat plate leading edge and forms a Blassius velocity 
profile where the velocity fluctuations are negligible. 
Further downstream the flow undergoes transition to turbulent 
flow [1] where the velocity inside the boundary layer 
exhibits irregular velocity fluctuations. The value of the 
Reynolds number, Re x , can provide an estimate of the laminar, 
transition and turbulent regions. The skin 


friction 
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coefficient for the laminar region of a boundary layer is 
smaller than that for the turbulent region. The most popular 
approach in wall friction drag reduction therefore has been 
to.__atjtempt_.to_. delay — the transition ^from_ laminar to turbulent 
flow [2] and maintain the laminar skin friction coefficient 
for as long as possible. Most practical flows are, however, 
turbulent and at large Reynolds numbers all flows will 
eventually undergo transition to turbulence. The main object 
of the following investigation is to reduce the skin friction 
drag of a turbulent boundary layer by altering its structure 
[3]. 

Visual investigation of turbulent boundary layers have 
shown them to contain at least three distinct types of eddy 
structure (figure 1) : 

i. Large eddies which dominate the outer region of the 
boundary layer and have an average length of 
approximately 1 . 6 <5^ [ 4 3 

ii. Typical eddies having a typical length of 200 wall 
units, y + , are encountered throughout the turbulent 
boundary layer [ 5 ]. 

iii. Well organized motions in the laminar sublayer can 
form low speed streaks [6]. 

The large scale eddies have an average axial length 
of about 1.6^ [4] and cause a peak velocity fluctuation of 
the order of 0.2U o where U 0 is the external velocity. The 
large scale structures are convected downstream with an 
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Figure 1. Typical flat plate boundary layer 
(flow visualization with smoke). 
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average velocity U C =0.93U 0 [7], and appear serai-periodically 

with an average frequency of appearance f = . [8], 

2 . 5 o 

These large scale structures have a limited life span. They 
-decay -after ; traveling- a "distance o-f- about" IOcT and lose - their 
coherence after about 4(5" [9]. Large eddies are associated 
with the production of turbulent energy and the production of 
Reynolds stress -pu’ v* , although the exact mechanism linking 
them to the sublayer skin friction producing phenomena is 
still unclear [4,6,9]. The reason for this is the lack of 
comprehensive experimental data in ~ the laminar sublayer. 
Consequently, " manipulation ‘ ~of the turbulent boundary layer 
structure could lead to reduced skin friction as described in 
the following sections. 


1.2. HISTORICAL REVIEW 


Early boundary layer manipulation [10] was performed by 
positioning mesh screens of height approximately equal to cT 
inside the boundary layer. This resulted in the destruction 
of the large eddies and a skin friction drag reduction of 50 
percent extending a distance of 100(fT to 150<fT downstream of 
the mesh screens. The drag of the device itself was 
excessive. Pursuing the same goal of eliminating the large 
eddies, Hefner et al [11] mounted horizontal and vertical 
plates of various shapes and sizes inside the boundary layer. 
These plates or elements interacted with and modified the 
large eddies, interrupting the turbulence production. The 
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reduction of the wall skin friction was attributed partially 
to the alteration of the outer turbulence scale. An average 
of 24 percent drag reduction for a distance of 45 device 
.heights was measured but no net drag reduction was observed. 
A study utilizing four parallel plate manipulators embedded 
inside the boundary layer was performed by Corks et al [12]. 
This method was effective in removing the large scale 
structures and a 30 percent skin friction drag reduction was 
documented at a distance of 70(5* downstream. Flow 
visualization recorded the breaking up of the large eddies by 
-the' ' parallel : -plates which resulted in a- thinner manipulated 
boundary layer. The streamwise turbulence intensity close to 
the wall of the manipulated boundary layer was found to be 
reduced. No net drag reduction was measured. In order to 
improve upon the net drag reduction Corks et al [13] mounted 
a pair of tandem flat plates within the boundary layer. They 
measured a 20 percent net drag reduction at some locations. 
The local skin friction drag coefficient did not relax to 
that of the undisturbed flow until a distance of 6d&. They 
also identified the possible mechanisms leading to such 
results as the suppression of the normal velocity component 
of the large eddies by the manipulators and the 

redistribution of the turbulent kinetic energy by the blade 
wake. Bertalrud et al [14] used thin transverse ribbons in 
tandem for various flow configurations. A 40 percent skin 
friction drag reduction was measured and this phenomenon 
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persisted for about 80 downstream of the manipulator. The 
net drag reduction was moderate. 

These experiments adopted similar boundary layer 
“manipulating ‘schemes, but the specific dimensions of chord 
length, height and tandem spacing varied considerably. In a 
review paper Hefner et al [15] presented the parameters for 
- an optimum plate -manipulation of a turbulent boundary layer. 
They stated that the flat elements should be thin and mounted 
as tension members. They estimated that the element chord 
length should be approximately equal to the boundary layer 
thickness,^ , and the outer element should not exceed a 
height of O.S^above the wall surface. It was also noted that 
the number of horizontal elements should not exceed two and 
that the plates operate most effectively in the range of 

momentum thickness Reynolds number from 3,000 to 6,000. In 

recent investigations [16], 10 to 30 percent net drag 

reduction was achieved with optimized flat plate manipulators 
in tandem under zero, mildly adverse and favorable pressure 
gradient flows. Most of the net drag reduction occurred 
within 50 boundary layer thicknesses from the manipulator. 

Airfoil shaped [17] large eddy break-up devices (LEBU's) 
were tested in a range of momentum thickness Reynolds numbers 
at the manipulator near 7,400. The results indicated that 
only symmetric airfoils in tandem were effective in 

manipulating the turbulence boundary layer and a 7 percent 
net drag reduction was recorded. Finally it should be 
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mentioned that the repeatability of net drag reduction 
experiments was not found satisfactory [15] and the 
microgeometry of the plate manipulators is critical with 
respect" to their effectiveness [18]. 

1.3. OBJECTIVE OF STUDY 

Large eddy break up devices (LEBU J s) constitute a new, 

promising method of obtaining drag reduction in a turbulent 
boundary layer. The destruction of the large scale eddies and 
the associated reduction in momentum transfer to the wall can 
result in lower skin friction drag. One of the suggested 
mechanisms, based upon the work of Liss and Usol'tsev 
[13*19>20], leading to large eddy elimination is 
schematically represented in figure 2. Large eddies 

(discrete vortices) impinging upon a LEBU plate generate 
unsteady lift forces upon it due to the effective angle of 
attack variation. Because of the unsteady circulation around 
the device, vortices are shed from its trailing edge. The 
oncoming large eddies in the boundary layer interact with the 
shed vorticity and are partially canceled. The focus of this 
dissertation is to enhance the LEBU effectiveness by exciting 
its trailing edge with acoustic waves phase locked to the 
large scale structure and thus influence the momentum 
transfer to the wall. 

An initial estimate of the required sound pressure level 
for an effective acoustic pulse was obtained by considering 
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Figure 2. Large eddy brake-up process. 
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the magnitude of the pressure perturbations at the near wake 
of a thin plate in inviscid flow. A sound wave of 
approximately 100 dB (Appendix 1) pressure level, arriving 
a-t -t-he LEBU -trailing : adg“e- together- ^ith the large eddies was 
considered adequate to influence the already manipulated 
boundary layer. 

Detailed skin friction measurements were obtained in the 
flow region downstream of a LEBU excited with acoustic waves. 
The data are compared with skin friction measurements of a 
simply manipulated flow, without acoustic excitation and with 
a plain flow configuration. The properties and the scales of 
motion in the flow regime downstream of the acoustically 
excited LEBU are studied. A parametric study based upon the 
characteristics of the acoustic input was pursued in addition 
to careful mapping of the drag reduction phenomenon within 
the acoustically manipulated boundary layer. This study of 
boundary layer manipulation has lead to improved skin 
friction drag reduction and further understanding of the 
turbulent boundary layers. 
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2. EXPERIMENTAL PROCEDURE 

2.1. BASIC EXPERIMENTAL APPARATUS 

— — The -experimental ~study„was:--conducted in the recently 

constructed NCSU Boundary Layer Wind Tunnel [21]. A general 
view of the wind tunnel is shown in figure 3- This facility 
features: a long test section (~7.0 m) to achieve high 

Reynolds numbers with low speeds, a velocity range of 4 to 
33 m/s and extremely uniform mean flow with turbulence 
intensities of less than 0.25 percent. A single LEBU plate 
zvas :: carefully fabricated from -stainless steel shimstock and 
mounted at a distance of 2.0 m from the leading edge of the 
test section, in the boundary layer (figure 4). A mounting 
:system for the LEBU was designed and built. The mount has 
accurate control over the LEBU's height above the test 
surface, the LEBU's angle of attack and tension. The LEBU 
extends through openings in the wind tunnel side walls across 
the full span of the test section. Each side of the LEBU 
mount is independently adjustable to assure precise alignment 
and eliminate possible vibrations of the stretched plate. 

The experimental arrangement to acoustically excite the 
LEBU trailing edge is shown in figure 5. A hot film probe is 
positioned upstream of the LEBU at the LEBU height. The 
signal from the sensor is obtained from an anemometer and fed 
into a specially designed processor whose electronic diagram 
is drawn in Appendix 3. The processor responds only to large 




Figure 3. General view of the NCSU Boundary Layer Wind 
Tunnel . 


LEBU BOUNDARY 
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Figure 4. Wind tunnel test surface 




Figure 5. Schematic represents tion of the acoustic 
excitation mechanism. 
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excursions of the hot film signal. The threshold voltage 
level, V r , to which the processor responds is adjustable. This 
adjustment allows selection of a lower limit of large scale 
- eddies - to' - which ; -the- processor - will — respond . Eddies with 
velocity fluctuations above this threshold level (reference 
voltage) are ignored while eddies with associated velocity 
excursions below the threshold level produce a response from 
the processor. The response is a signal with character and 
duration similar to the anemometer’s large fluctuations as 
shown on the dual trace oscilloscope display of figure 6. The 
response of ' the processor must also be' time delayed and 
amplified. The time delay, t, is set equal to the large eddy 
convection time [7] between the upstream hot film sensor and 
the LEBU trailing edge. - The delayed processor response is 
then amplified to produce an acoustic input capable of 
influencing the LEBU wake. 


The processor response is then used to drive acoustic 


waves from a 

location on 

. the floor 

of 

the wind tunnel 

below 

the trailing 

edge 

of 

the LEBU. 

The 

acoustic waves 

can 

therefore be 

made 

to 

arrive at 

the 

LEBU trailing 

edge 


together with the incident large eddies. In this manner the 
vorticity shed from the LEBU, which apparently influences the 
eddy cancellation, can be modified. The pressure pulse device 
consists of an 8 inch woofer speaker which is seal mounted on 
the wide side of a conical channel. The speaker apparatus is 
housed in a lined box to restrain noise signals from entering 
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the test surface by spurious paths. 

2.2. OPERATING CONDITIONS 

-i-fChe”-— wind- tunnel— is adjusted to- produce, a zero pressure 

gradient flow along its length. The variation of the free 

stream turbulence intensity in the test section is shown . in 

figure 7. For the operating velocity range a 0.25 percent 

turbulence intensity was measured. The boundary layer is 

forced turbulent by a coarse sandpaper trip mounted 3 cm 

downstream of the test surface leading edge. The length of 

- -the- - -boundary layer -trip -is 21 xra -whrcti -quickly establishes a 

fully developed turbulent boundary layer with a complete 

range of length scales as can be seen from the photograph in 

figure 1. The wind -tunnel unit 'Reynolds number, -^2, was set 
6 —1 

at 0.7x10 m. The corresponding typical mean flow velocity was 

approximately 11.0 m/s and the Reynold's number based upon 

the momentum thickness , Q , at the LEBU location was about 

3,100. The boundary layer thickness, (5" , at LEBU location was 

estimated as 4.1 cm. This was determined using the definition 

of the boundary layer thickness as the height at which 

-^-=•995, where u is the local velocity and U Q is the 
Uo 

velocity of the external flow. The LEBU chord length was 
-0.90(5*, its height above the test surface was -0.80(5" and its 
thickness was 0.005(5" according to the recommendations of 
reference [ 15] . 

This LEBU configuration is an improved version of 
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several manipulators that were tested. It is described in 
detail because a major portion of this study was spent upon 
it. This plate configuration was the fifth tested (LEBU #5) 
and=rwas the most successful. The “configurations previously 

tested also used a wind tunnel unit Reynolds number of 0.7x 

6 -1 

10 m , were located at a distance of 2.0 m from the test 
section leading edge and .are -described in section 3.6. 

The upstream eddy detector hot film probe was positioned 

at the midspan of the test surface at a height of 0.855" and 

at a distance, s, of 2.455" upstream the LEBU trailing edge as 

shawn- ~in -figure 5. This probe, used to detect the incident 

large eddies, was positioned as close as possible to the LEBU 

trailing edge because the large eddies lose their coherence 

after a distance of about 45"f4]. The location of the large 

eddy detector probe must be, however, far enough upstream in 

order to prevent a feed back from the acoustic pulses. The 

axis of the probe was inclined at 45 degrees with respect to 

the test section vertical plane of symmetry (figure 4) in 

order to avoid probe interference with the downstream flow 

field. The processor delay time was set at t= — - ; a 

0.93 u 0 

typical value of t was 10 msec. 

The pressure pulse mechanism was mounted under the test 
surface at midspan below the LEBU device trailing edge. The 
acoustic input was directed to the LEBU trailing edge through 
a circular opening (3/8 inch diameter) covered with fine 
screen flush mounted to the test surface. 
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A first approximation of the required sound pressure 
level of the acoustic waves was obtained by considering the 
boundary layer pressure perturbations near the trailing edge 
-rof a~flat plate .-According- to - detailed calculations (Appendix 
1) a typical value of about 100 decibel (dB) pressure 
perturbation was predicted at the LEBU trailing edge. This is 
not the pressure fluctuation generated at the flat plate 
trailing edge when it encounters an oncoming average large 
eddy, but rather the typical pressure fluctuation in that 
location for a flat plate in inviscid free stream flow. It 
-was- assumed that sound pressure levels of this magnitude 
would be necessary to have some effect upon the trailing edge 
flow. A power amplifier was utilized to deliver the required 
power to the pulse mechanism. 

The pressure pulse mechanism was bench tested with the 

6 —1 

wind tunnel operating at a unit Reynolds number of 0.7x10 m 
using the output of the upstream hot film probe. The speaker 
mechanism was mounted outside of the wind tunnel under a 
wooden baffle plate with an acoustic wave port identical to 
that in the wind tunnel test section. The time averaged Root 
Mean Square (RMS) power spectra of the acoustic output at the 
LEBU's height was measured in a hard wall room under zero 
flow conditions. Figure 8 shows the RMS broadband spectra 
from 0 Hertz to 2,000 Hertz of the acoustic pulse emanating 
from the floor pressure port. The acoustic power is 
concentrated over a narrow region from 60 Hertz to 200 Hertz 



RMS 96.4 dB 



Figure 8. Narrow band acoustic signal of the pressure 
pulse at LEBU's height. 
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approximately at a level of 100 decibels. The sound pressure 
level at the LEBU’s height integrated from 0 Hertz to 2,000 
Hertz was 96.4 decibels. The large eddy passing frequency was 
-expec-ted -to-be - approximately 1 tQJiectz , hence the bulk of the 
acoustic power occurred at the expected frequencies. The 
remaining part of the spectra consists of the harmonics of 
the acoustic pulses, the background noise, "ground dip" 
phenomenon and the harmonics of the wind tunnel fan blade 
passing tone. A 1/2 inch Bruel and Kjaer condenser 
microphone was utilized for these measurements. 

2.3. DETERMINATION OF THE SKIN FRICTION COEFFICIENT , c f 

Measurements of the boundary layer profile downstream of 
the manipulator plate were obtained using a small boundary 
layer Pitot probe with an outside diameter of 0.5mm. A 
remotely controlled traversing mechanism was used to assist 
with the precision probe positioning. The traversing device 
provides probe movement over the Y-Z plane in the wind tunnel 
test section (figure 4) with a positioning accuracy of +0.1 
mm. The pressures from the boundary layer probe were measured 
with a Validyne differential pressure transducer. The static 
pressure was obtained from static taps in the test surface. 
Approximately 200 samples were obtained and averaged at each 
point and the velocity profile was determined by 40 such 
points distributed over the thickness of the boundary layer. 
A Digital 15 bit resolution data acquition system in 
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conjunction with a Digital Computer were utilized in order to 
process the measured data and monitor the wind tunnel' s 
operation. Variations in mean velocity of more than 1 percent 
and.- mean- flow temperature of more than 2 percent 
automatically stopped the acquisition process, so that the 
exact conditions could be reset. 

The momentum thickness , Q , is calculated from its 
definition (list of symbols) by numerically integrating the 
velocity profile over the boundary layer thickness. For zero 
pressure- gradient two-dimensional flow the Von Karman 
^Momentum •-Integral -Equation reduces to: 



( 2 . 1 ) 


where, c f is the local skin friction coefficient and x is the 
axial length. The momentum thickne'ss, Q, is determined as 
described above for a plurality of positions downstream of 
the LEBU. The values can be least square fitted to a power 
curve of the form: 


0=ax b 


( 2 . 2 ) 


Hence, the skin friction coefficient can be estimated by : 
c f =2abx b_1 (2.3) 

An alternative calculation of the skin friction coefficient 
can be obtained by applying a central finite difference 
scheme upon equation 2.1 

The previously described indirect method of measuring 
Cf , although time consuming, has certain advantages compared 
to other methods of wall shear stress measurement [22]. Skin 



friction gauges have problems with the nessecary gaps and 
possible misaligneraent related to the floating element. 
Preston tubes, Stanton gauges and sublayer fences cause flow 
obstruction, and measurements are heavily dependent ^upon 
geometry and positioning. The method used in this study was 
- the --most -economically sound and the most directly related to 
proven theory. 

2.4. VERIFICATION OF THE BASIC FLOW CONFIGURATION 

The .momentum thickness ,0, of the turbulent boundary 
... layer. ... is of primary interest .in -the present investigation 
because it is used for determining the skin friction 
coefficient. Figure 9 shows measured and theoretically 
predicted -values of Q for the -plain flow case versus 
downstream distance x. The solid line represents values of 
momentum thickness predicted by Head's method [23]. The 
experimental values of Q were obtained by averaging 3 values 
across the span of the test surface. Head's method was 
applied for zero pressure gradient flow and initialized with 
a momentum thickness predicted empirically at the end of the 
sandpaper trip (Appendix 2.1). The initial shape parameter, 
H=1.375, was selected as the averaged measured value. The 
shape parameter H for a boundary layer is defined as the 
ratio of displacement to momentum thickness. Measured values 
of H at different axial locations are listed in table 1 along 
with the average. The initial or starting values in Head's 
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Table 1 . Measured shape parameter, H, versus 
distance, x (-^=0.7x10 6 m _1 ). 

x [ra] H 

2.00 1.379 

2.25 1.332 

2.50 1.435 

2.75 1.375 

3.00 1.361 

3.25 1.363 

3.50 1.366 

3.75 1.343 

4.00 1.356 

4.25 1.353 

4.50 1.357 


axial 


average H=1 . 375 
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method are relatively unimportant except for the first few 
points of the calculation. Agreement between the predicted 
and measured values of appeared satisfactory with a 
-maximum-deviation of ~3 -percent; — Thor momentum thickness , Q , 
at 9 spanwise locations across the test section centerline at 
2.0 m and 3*5 m from the test surface leading edge is shown 
-on- figure -10. The maximum deviation from the mean was 4 per- 
cent. 


„ -0-2 6 8 -0.6 78 H 

C — 0.2 4 6 R eg 1 0 


The local skin friction coefficient for a plain flow 
configuration was determined according to the procedure 
suggested by Bradshaw [24]. This method explained in Appendix 
2.4 determines c^ from a measured boundary layer velocity 
profile. These measured results are compared with the axial 
skin friction coefficients predicted by the empirical formula 
of Ludwieg and Tillmann [1], 

(2.4) 

in figure 9. The shape parameter H and momentum thickness Re 
in the above equation were predicted with Head's method. A 
maximum deviation of 2 percent was observed between the c f 
predicted by Ludwieg and Tillmann' s formula and the Cf 
determined by Bradshaw's method. Although Bradshaw's method 
predicted c^ very well for the plain flow, it was not used in 
other flow configurations because it is based upon the 
previous knowledge of the logarithmic law ; the form of the 
logarithmic law for a LEBU manipulated boundary layer has not 
been established yet. 






The measured velocity profile for the plain flow 

configuration was compared to a theoretically predicted 

profile using a modified Coles' boundary layer (Appendix 

-2.2). The -velocity profile was calculated for an axial 

location 2.5 m from the test section leading edge. The free 

stream velocity was set at 11.35 m/s to satisfy the unit 

0 

Reynolds number requirement of 0.7x10 m and the resulting 
theoretical and experimental profiles are shown in figure 11. 
The maximum deviation between predicted and measured velocity 
values was 2 percent. 

— Measurements of -the fluctuating -axial velocity component 
were obtained across the turbulent boundary layer at an axial 
location of 2.0 m with a hot film probe. These data ar 
-plotted in figure 12 -and compared with classical measurements 
of Klebanoff [1], The trend of the data agrees fairly well 
although the magnitude of the measured data is lower than 
Klebanoff ' s classic measurements. Such deviations should be 
expected since Klebanoff' s data are for a Reynolds number of 
10 7 whereas the measurements of figure 12 were obtained with 
a Reynolds number of 1.4x10 6 . A method to scale the measured 
data to correct for the difference in Reynolds numbers was 
not available. 

The boundary layer velocity profile for a LEBU 
conf iguration was predicted by superimposing the turbulent 
wake velocity distribution from a thin plate upon a plain 
flow profile. The detailed calculation procedure is described 




Predicted and measured mean velocity profil 
for the plain flow. 




Figure 12. Measurements of the streamwise velocity 
across the boundary layer. 
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in Appendix 2.3 and the predicted and measured profiles are 

shown in figure 13. Their agreement was found satisfactory and 
a 4 percent maximum difference was measured. 

The structure of the turbulent boundary layer was 
investigated by using flow visualization with the Smoke-Wire 
method of Raspet and Moore [25]. The smoke was produced by 
- vaporizing - oil from a- fine stainless steel wire by impulsive 
resistive heating. The smoke wire was mounted vertically in 
the test section with the tunnel operating at low speed. 
Large scale eddies shaping the boundary layer intermitancy 
region and typical eddies throughout the layer can be 
observed in figure 1. 

The mean free stream velocity was uniform in the 
spanwise direction within 0.2 percent [21] and momentum 
thickness across the test section span varied within 4 per- 
cent of the mean value. These small deviations combined with 
the satisfactory agreement found between predicted and 
measured flow characteristics allow the assumption of a two 
dimensional ergodic turbulent boundary layer with the large 
eddy break-up device removed. When the LEBU is installed the 
basic velocity profile downstream of the device is also as 
expected. The wake of the plate appears normal and the slight 
increase in the velocity profile, at approximately y=1.0 cm 
in figure 13 » may be attributed to the viscous drag reduction 
capability of the LEBU which manifests itself as a reduction 
in momentum thickness. 




LEBU ff 2 configuration 
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3. RESULTS AND DISCUSSION 

The experiments with acoustic excitation of the LEBU 
were performed in order to: 

i. Determine if proper acoustic excitation can enhance 
the effectiveness of a boundary layer plate 
manipulator . 

" ii .'" Investigate "the ' "properties ~of the resulting flow 
field downstream of. the LEBU. 

iii. Search for the optimum operating parameters of the 
experiment. 

iv. Examine the effectiveness of the acoustic excitation 
with different experimental conditions. 

■ _1 3.1. BASIC RESULTS 

3.1.1. SKIN FRICTION DRAG REDUCTION 

The effect of the acoustic pulse upon the plain flow 
was first examined and the momentum thickness growth for the 
unexcited and the excited plain flows is compared on figure 
14. The acoustic input caused no significant alterations in 
the flow and the average value of 0 for the acoustically 
excited configuration was found to be higher than the 
unexcited configuration by only 0.5 percent. This slight 
increase in 0 is due to an expected small increase in mixing 
near the acoustic source. Figure 14 shows that the acoustic 
input has no significant effect on the flow when the LEBU is 



Unexcited 



Figure 14. Momentum thickness versus axial distance 
,X» for the excited and unexcited plain 
flows. 
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not installed. 

When the thin plate was installed turbulent boundary 
layer velocity profiles were obtained at raidspan locations 
downstream of - the .LEBU manipulator . The momentum thickness , 9 , 
was calculated as described in section 2.3 and plotted versus 
downstream distance ,x, from the test surface leading edge. 
Figure 15 shows the axial variation of for the following 
configurations : 

i. Plain flow with the upstream sensor installed 
ii* The LEBU configuration with the upstream sensor 

installed but not operating. 

iii. The LEBU conf iguration acoustically excited. 

The momentum thickness of the manipulated flow is 
greater than the plain flow case near the LEBU because of 
momentum loss imposed by the embedded thin plate. The 
momentum thickness near the LEBU trailing edge ( 4 boundary 
layer thicknesses downstream) is within 5 percent of the 
sum of the momentum thickness of the plain flow case and the 
momentum thicknesses of the laminar boundary layers that 
develop on the upper and lower surfaces of the plate. 

The acoustically excited case produces an even greater 
momentum loss at this location. If the acoustic input was 
merely tripping the laminar boundary layer, one might hope to 
repeat the above superposition using a calculation for 
turbulent boundary layers on the LEBU plate. The measured 
increase in Q with the acoustic input is, however, 
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Figure 15. Momentum thickness, o, versus axial distance 
,X, for various flow configurations. 
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approximately ten times the expected increase that could be 
attributed to turbulent boundary layers. These results 
suggest that an additional mechanism may be active during the 
acoustic excitation; a mechanism _which alters the trailing 
edge flow from the blade in a manner other than merely 
tripping a laminar layer to a turbulent boundary layer. 

Further downstream at x = 2.6 ra (15 S' from the LEBU) the 
momentum thickness of the acoustically excited configuration 
relaxes to that of the unexcited LEBU but grows at a slower 
rate. -The same trend is -noted at x=3.S m where the momentum 
irfehiakness^xof._the -acoustical ly^excited:f low becomes less than 
the plain flow conf iguration . This a clear indication of 
reduced skin friction coefficient and net drag reduction. 

---The-- -approximate ^variation of the skin friction 
coefficient, c f , with the axial distance from the test surface 
leading edge is shown in figure 16. All measurements were 
obtained at the test surface mid span. This location was 
aligned with the upstream sensor and the acoustic wave port. 
The numerical values of the skin friction coefficient were 
calculated according to the procedure described in section 
2.3. All least square curve fits had a correlation 

coefficient above 0.99. The LEBU conf iguration and the 
acoustically excited flow were compared to the plain flow 
(hot film sensor installed). According to the numerical 
values of the skin friction coefficient, acoustic excitation 
can enhance the LEBU's effect by reducing c f between 


PLAIN FLOW 



Figure 16. Variation of the skin friction coefficient 
t C, , for each conf iguration . 
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approximately 7 and 13 percent with respect to the unexcited 
LEBU configuration. Comparison with the plain flow shows a 
reduction of the wall shear stress between 9 and 15 percent. 

3 .T. 2. -"SXITr TRICTION DRAG REDUCTION MECHANISM 

Space-time cross correlation functions (CCF) ■ were 
obtained using the hot film eddy detector probe upstream of 
the LEBU plate as a fixed reference location and another, 
probe at various downstream locations as indicated in figure 
17» A dual channel anemometer was utilized and the hot film 
p robe — signals .were processed^: j*ith : :a- NICOLET SSOB FFT 
analyzer. Figure 13 shows peak values of the space-time cross 
correlation functions versus probe separation distance,^, at 
various boundary -layer -heights .--far away from the LEBU, in 
the downstream direction, the residual acoustic far field 
produces higher cross correlation peaks whan the acoustic 
excitation is applied. When the moving probe is located near 
the LEBU's leading edge the excitation causes significant 
increases in the cross correlation functions because the 
probes are exposed to both, the large scale eddies and the 
acoustic field. At the trailing edge of the plate in the 
outer region of the boundary layer the flow without 
excitation maintains a maximum CCF approximately equal to 
that of the leading edge. The wake of the LEBU in this case 
is clearly correlated to the incident flow. When acoustic 
excitation is applied the downstream flow is correlated to 



moving 

PROBE 



Figure 17. Schematic for velocity cross correlation 
measurements across the LEBU. 
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tha upstraam point significantly less. The reduction is more 

y 

obvious at the outer region of the boundary layar ( 4r= 0.95 ) 

0 

where tha large eddies dominate. This indicates that acoustic 
excitation enhances the destruction of the large eddies which 
are sensed at the upstream location. It is assumed that the 
impinging sound wave at the LEBU trailing edge sheds a 
vortax of opposing rotation to the oncoming large eddies 
resulting in enhanced large, eddy cancellation . 

At downstream locations deeper inside the boundary layer 

y 

-j=r =0 . 60 ) the reduction of the CCF, when acoustics 
O 

"are applied, is moderate because -large- eddies appear mainly in 
tha outer part of the boundary layer. Tha reduction of the 
CCF when the moving probe is at the LEBU's height signifies 
- that -the acoustic pulsas assist- the -manipulator in canceling 
large eddies. All time delay values associated with the cross 
correlation function peaks correspond to expected large eddy 
convection times and not to the acoustic wava propagation 
time . 

The turbulence intensity of the low frequency part of 
the turbulence spectrum was examined downstream of the plate 
manipulator. A hot film probe was positioned slighty above 
(1.5 mm) tha LEBU's height and the turbulence spectrum from 0 
to 2,000 Hz was measured. Figure 19 compares tha variation of 
the root mean square (RMS) of the fluctuating axial velocity 
component u’ for the acoustically excited flow to the 
unexcitad LEBU conf iguration at various downstream locations. 


( = 0 . 80 , 
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Figure 19. Relative variation of the low frequency RMS 
turbulence close to the LEBU trailing edge. 
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The turbulence level over 2 KHz is considerably reduced when 
acoustics are applied and the effect persists for a distance 
of 8 (^ downstream of the LEBU. When the low frequency (0 to 
- 200 — Hz - ) “part of the signal is examined the level of the 
turbulence intensity appears further reduced. This range of 
frequencies approximately corresponds to the large eddy 
passing frequency which is estimated to occur near 110 Hertz 
and dominate the low frequency spectrum. The relatively 

higher frequency (200 to 2,000 Hz) part of the turbulence 

spectrum remains largely unaffected by the acoustic pulses. 
This indicates that only the large wavelength motions of the 
turbulent boundary layer are affected when the flow is 
acoustically excited. The measurements of the RMS velocity 
fluctuation, as for the previous CCF data, were limited to a 
distance of about 10^ downstream the LEBU. 

To better understand the influence of the acoustic 
excitation on the flow manipulator, a pulsed smoke wire was 
used to obtain flow visualization photographs. The smoke wire 
was mounted downstream of the manipulator at the midspan of 

the tunnel. The fine wires ran vertically from the wind 

tunnel floor through the boundary layer. Figure 20 shows 
typical results for the three cases of concern. The large 
eddy break-up device reduces the large scale structure in the 
boundary layer and acoustic excitation of the device causes 
further reduction in coherent motion. The results were 
obtained at approximately 35(^ downstream of the manipulator. 



ORIGINAL PAGE IS 
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Figure 20. 

I 


Boundary layar flow visualization for tha 
(A) plain case (B) LS3U manipulated and (C) 
acoustically excited LS3U configuration. 
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The sound waves at the LEBU trailing edge apparently shed 
vorticity which helps cancel the oncoming large eddies. 

The possibility of vertical oscillation of the thin 
ribbon when "hit" by acoustic waves was examined. The goal 
of this study is to control the large eddy break-up process 
without moving the LEBU. M_ot.ion_of._the device could generate 
more disturbances and altar the boundary layer on the device, 
thus :_altering - the -momentum thickness downstream. A small 
accelerometer was mounted on the LEBU above the acoustic wave 
port. Vertical displacement measurements were obtained with 
no flow; the pressure pulse device was excited with a 
previously recorded anemometer signal from the turbulent 
boundary layer. The measured vertical displacement when 
acoustics were applied could not be distinguished from the 
background electronic noise. The reduced correlation of the 
flow across the LEBU and reduction of RMS velocity 
fluctuations with acoustics, can be only attributed to the 
interaction of the acoustic input and the fixed LEBU, not the 
motion of the manipulator device. This is evident because the 
acoustic pulses had no effect upon the the plain flow; also 
they cause no vertical motion of the plate at zero flow 
conditions. 

According to these results an acoustic pulse of the 
proper sound pressure level and phase locked to the convected 
large scale structures can enhance the effectiveness of a 
plate manipulator. Large scale eddies contribute 
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approximately 50 percent of the turbulent energy and 30 per- 
cent of the Reynolds stress when the outer part of the 
boundary layer is considered [9]. Turbulent "bursting" in the 
sublayer accounts for most of the turbulence and Reynolds 
stress close to the wall [26,27]. If the passage of the large 
eddies is related to the triggering of "bursting" events in 
the sublayer, then their enhanced destruction can lead to 
further reduction of skin friction drag downstream of the 
LEBU. Energy exchange between the mean flow and turbulence is 
governed by the _ dynamics of, the large eddies [23] which 
--extraci- kinetic. energy from -the- .mean flow. The break-up of 
the energy cascade from the very large eddies to the small 
energy dissipating structures can lead to reduced momentum 
_ transfer- to_ the- wall- .and .reduced skin friction drag. 

3.2. SPREADING EFFECT 

It is clear from figure 15 that increased drag reduction 
due to acoustic excitation persists downstream for a 

considerable distance. Measurements of Q across the span of 
the test section indicated that the beneficial effects of 
excitation spread slowly in the spanwise direction as the 
flow convected downstream. 

3.2.1. SPANWISE MAPPING DOWNSTREAM OF THE LEBU 

The spreading of the acoustic excitation effect 
downstream of the LEBU was mapped indirectly. It has been 
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observed (section 3.1.2) that the output signals of two hot 

film probes across the LEBU at midspan of the test surface 

are less correlated when acoustic excitation is applied. The 

cross correlation-function~(-SCF^ between the upstream fixed 

sensor and another positioned at various spanwise locations 

downstream of the blade was measured. The ratio r r 

c CCRL) 

was calculated; LA signifies the acoustically excited LEBU 
configuration and L signifies the LEBU configuration without 
excitation. Figure 21 shows a schematic top view of the test 
surface downstream of the LEBU with expanded z coordinates. 
The solid symbols indicate points on the test surface where 
the ratio r c is smaller than unity because the acoustic waves 
enhance the large eddy destruction. The open symbols indicate 
points where the ratio r c is greater than unity because the 
acoustic pulses make the probe signals more correlated. These 
measurements could not be extended further downstream because 
the large eddies have a limited life span and decay after 
traveling a distance of approximately 10(5*. A slight 
spreading of the acoustic excitation effect at a half angle 
of about 2 degrees can be seen in figure 21. This spreading 

is verified by spanwise measurements of momentum thickness. 

L A) 

Figure 22 shows the variation of the ratio Tq ~ 

with the spanwise distance, z, from the test surface 
centerline for various axial locations. The variation of 
momentum thickness , Q , (section 3.1.1) measured at the test 
surface centerline needs to be reconsidered. According to 
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Figure 21. Spanwise mapping of the peak cross corre- 
lation downstream of the LEBU (solid symbol 
indicate less than 1.0). 
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Figure 22. Spanwise variation of momentum thickness ,0 
for the excited and unexcited LEBU at 
three axial locations. 
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these data the ratio rg is greater than unity at axial 
locations close to the LEB'J trailing edge. Further downstream 
at a distance of 6<f to 8 ( 5 " from the LEBU trailing edge the 
rail 0 -rg becomes smaller -than unity, — remaining as such for the 
rest of the test section. Returning to figure 22, at x=2.lS 
m, the spanwise measurements of Q indicate a narrow region 
where is greater than unity. Further downstream, at axial 
distances of 22(5" and 47(5" from the plate manipulator, the 
ratio rg was found to be smaller than unity over a broader 
extent in the z direction. These data show a tendency for the 
drag - reduction - phenomena to spread spanwise at a half angle 
of about 2 degrees. The experiment was limited to a distance 
of 50 boundary layer thicknesses from the LEBU because most 
of the previous investigations indicated that a manipulated 
turbulent boundary layer relaxes to the undisturbed flow 
after 50 to IQO^. Furthermore, the length of the test section 
did not allow good measurements beyond this point. The small 
value of the spreading angle is attributed to the mixing 
process of the surrounding turbulent boundary layer. This 
angle also agrees with the results observed for the spanwise 
spreading of the reduction in the space-time cross 
correlation near the LEBU. 

The turbulence intensity, T, of the axial flow 
velocity component, with and without excitation was measured 
throughout the boundary layer at three axial locations. 
Figure 25 shows the variation of the parameter r T = - ^ L . A ^ , 
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Vu y 

where T=— = — , versus the boundary layer height — at a 

spanwise location 2cm off the test surface centerline, z=+2.Q 

cm. 

. At— x=2. OB_ra^_ _ outside of_ the.-. narrow "wedge" of acoustic 

influence, the average value of the parameter r T is 1.003. 
Acoustic waves in this case make the excited flow more 
perturbed. At x=2.S3 m, where this point is starting to be 
affected by the acoustic influence, the average ratio is 
1.001. Well inside the cone of acoustic influence, at an 
axial distance of 47 $ from the LEBU, the average value of 
-the . parameter . is~3mallarr:T;han unit.yV3 It should be noticed 
that the turbulence intensity level includes all frequencies, 
not just the lower frequency range which is primarily 
influenced by the - acoustic - input. The differences in 
turbulence intensity are small when acoustics is applied 
leading to a random oscillation of r T around unity. The 
reduction of turbulence intensity seems to exist however, 
inside the narrow "wedge" of acoustic influence. 

3.2.2. EXPLANATION OF THE SPREADING EFFECT 

A recently proposed model of turbulent boundary layer 
coherent motion [29] was examined and a qualitative 
explanation of the spreading acoustic influence was based 
upon it. According to this model, large scale motions are 
initiated by the instability of the main flow and rotate in 
the direction of the mean flow shear. The sublayer bursting 
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events are caused by the interaction of the semi-periodic 
large scale motions with the wall. The burst is bound by a 
pair of counter rotating vortices which are inclined forward 
and the resulting ejection and sweep are flows induced by 
these vortices. A schematic representation of the suggested 
model is shown on figure 24a. 

Using reference [33] applied .to the particular 
experimental set-up, an . approximate sketch of the resulting 
flow field from a single pressure pulse can be made as shown 
in figure 24b. This schematic is a side view of the flow 
fireld- -Induce' d-^by -an acoustic - pulse '-Impinging upon the LEBU as 
seen from a downstream location. The acoustic pulse is 
spread through the air with the vortex ring moving away 
from the acoustic input port. The resulting jet flow has a 
spanwise velocity component w' which is also augmented by the 
stagnation of the pulsating air upon the LEBU. The moving 
media causes a diverging motion of the generated vortices 
through the downstream boundary layer. 

A top view of the flow across the flat plate manipulator 
is shown on figure 25. Upstream of the LEBU the large scale 
structures appear according to the model at a random rate and 
spacing. Only the eddies that are convected on the centerline 
are detected and canceled with acoustic waves. These acoustic 
waves partially stagnate upon the flat plate and induce a 
vortex "street" which is spread into the boundary layer 
downstream. These vortices interact out of phase with some of 
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Figure 25. Top view of suggested acoustic spreading 
mechanism. 
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the large eddies which were not detected by the hot film 
probe but were partially canceled by the manipulator. This 
cancellation proceeds downstream within a narrow angle 

because the spanw_ise velocity c.omp.onent w' of the induced 

secondary flow is assumed far smaller than the axial velocity 
component u. Some of these vortices eventually decay because 
of strong mixing in the turbulent boundary layer. It should 
be noted that the angle of acoustic influence was reduced at 
higher free stream velocities which supports this proposed 
explanation of how the' spanwise spreading of the phenomenon 
may occur. 

An alternative explanation can be provided based upon 
the assumption that the large eddies and the sublayer 
phenomenaof "burst" and ' "sweep" are closely interrelated and 
preserve each other. It is assumed that a burst of fluid is 
counteracted by large scale motion in order to preserve 
continuity at a point of the flow field. Also, the convected 
large eddies, being limited spanwise, trigger bursting events 
at other nearby locations in the z direction. These events 
account for most of the skin friction production. 
Consequently, the break-up of large eddies will spread 
spanwise interrupting the chain of events that might link the 
coherent motion to the sublayer bursting and sweeping of 
fluid. 
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3.3. OPTIMIZATION OF THE EXISTING CONFIGURATION 


Since acoustic excitation reduces the correlation 
between two points across the LEBU, this was utilized to 
optimize the various parameters which control the function of 
the analog time delay, processor . _ These parameters are the 
voltage reference ,V r , the time delay, t, and the pressure 
level of the acoustic pulse. 

The downstream hot film sensor was positioned at a 
distance of 7cra ( 1 . 7 &) from the LEBU trailing edge, and at a 
height of 3-4 cm above the wall on the centerline of the test 
surface. The upstream hot film sensor was positioned at the 
usual distance of 2.5(5" from the LEBU trailing edge. The 
ratio r C = ' cc ' ^g ^ s plotted versus the relative processor 


reference voltage, 


v a 


where V r is the sat reference 


voltage and V a is the upstream anemometer voltage output. 

The values of r c ( open symbols) displayed in figure 25 

Vr 

indicate a minimum at — — =0.970. The turbulence intensity 

v a 

level at the upstream sensor' s height was 2 to 3 per cent. 
The optimum range of the reference voltage is vary narrow (10 
millivolts) and it was repeatedly noticed during the 
experiments that V r should be set just below the AC RMS level 
of the detector probe output, V a . 

The closed symbols of figure 25 represent the averaged 
number of processor pulses per second for different values 
of the ratio . The average processor response frequency 

v a 

was determined by counting the number of output pulses from 
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Figure 25 . Relative peak correlation , r c . and number of 
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the processor over an interval of 200 sec. For the local 
boundary layer thickness of about 4 cm and free stream 
velocity of about 10.9 m/sec, a value of approximately 110 
Hertz is predicted [9] for the large eddy passing frequency. 
This frequency is very close to that measured when V r is set 
to its optimum value. Thus optimization of V r is again 
associated with the large scale structure in the boundary 
layer. 

The acoustic excitation remained effective whan the 
reference voltage was ..set _be_low the optimum value. In this 
case only very large eddies were recognized and the number of 
acoustic pulses per second was reduced. When the relative 

Vr 

voltage was set at 0.91 _the .processor was unable to 

v a 

trigger the pulse mechanism. On the other hand, the acoustic 
excitation becomes lass effective when the reference voltage 
is increased above the optimum range. In this case the 
acoustic excitation is triggered by additional smaller 
eddies. The number of acoustic pulses per second was 
increased to almost double the predicted number of large 
eddies per second. Increased correlation measured between the 
two points across the LSBU is than caused by acoustic waves. 
The optimum value of the reference voltage, V r , depends upon 


the anemometer voltage 

output , 

V a , which 

varies with 

the 

tunnel free stream ve 

locity. 

Continuous 

monitoring 

and 

periodic adjustment 

was nece 

ssary to 

compensate 

for 


temperature variations of the flow. 


61 


The procsssor time “delay, t, required to allow for a 
detected large eddy to be convected to the LEBU trailing edge 
was also optimized based upon the reduction of peak 
correlation across the LEBU when acoustics is applied. On 
figure 27_ the ratio. r c = is plotted versus the 
processor time delay t. The downstream probe remained at the 

-same position and the reference voltage was set to the 

optimum value. According to the data, a minimum value of the 
ratio r c can be distinguished corresponding to a time delay of 
10 msec. For the particular__tunnel ..velocity setting of 10.5 
m/sec, the time delay, t=- ^ , was found to be 10 msec. This 
reinforces the concept of assisting the large eddy 
cancellation process with an acoustic pulse which arrives 
simultaneously at the trailing edge of the plate manipulator. 

Finally, the time averaged sound pressure level of the 
acoustic waves impinging upon the thin plate was optimized 
using the same phenomenon of minimized CCF when acoustic 
excitation is applied. Measurements of the time averaged 
pressure level of the acoustic wave were obtained at 
"simulated" experimental conditions. The acoustic pulse 
mechanism was mounted upon a wooden baffle plate that had an 
acoustic port identical to that of the test section. The 
pressure pulse mechanism was triggered by a prerecorded 
output signal from the upstream sensor at real flow 
conditions. Sound pressure level measurements were obtained 
at the LEBU' s height above the pressure port in a hard wall 
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Figure 27. Variation of the relative peak correlation, 
r c , with time delay, t. 
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room with no flow. The microphone was positioned at grazing 

incidence and the reference voltage, V r , and time delay ,t, 

were set at their optimum value. In this manner the sound 

pressure level at the LEBU was determined for various power 

amplifier settings. Measurements of the CCF were then taken 

inside the tunnel with real flow conditions. A hot film probe 

was positioned at midspan of the test section at a distance 

of 1.7(5" from the LEBU trailing edge and at a height of 0.82fi 

from the wall. CCF measurements were obtained comparing the 

.LEBU.__.and the acoustically excited flow configurations for 

different settings of the acoustic wave pressure level. The 
cc f(i 

ratio r r = ^ ) u x ’ versus the measured narrow band SPL at the 
c CCF(L) 

LEBU's height is plotted on figure 28. Once more a minimum 
value can be observed, which indicates that an acoustic pulse 
with a SPL of about 102 decibels is most effective for 
canceling the large eddies. 

The cross correlation across the LEBU remained reduced 
for a range of 2 dB around the optimum value but for values 
of the sound pressure level above 105 dB the acoustic wave 
field increases the correlation across the plate. The value 
of the sound pressure level required for optimum large eddy 
manipulation is of the same order of magnitude as the 
predicted and previously applied pressure pulse of 100 dB. In 
the beginning of the study it was assumed that a pressure 
pulse of the same order of magnitude as the pressure 
perturbations at the near wake of the thin plate (Appendix 1) 



no flow 



Figure 29. Variation of the relative peak correlation, 
r c , with time averaged SPL of pressure puls 
at LEBU' s height. 
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might benefit the large eddy break-up process. 

The observed phenomenon of reduced correlation at two 
points across the LEBU device whan the flow is acoustically 
excited can be utilized to optimize the operating parameters 
of the acoustic input. The cancellation of the large eddies 
is maximized when the reference voltage is set to produce 
approximately' the same number of acoustic pulses per unit 
time approximately as the number of predicted large eddies 
per unit time. The flow is least correlated ( enhanced large 
eddy cancellation)when the time delay is set equal to the 
large eddy convection time from the detector probe to the 
LEBU trailing edge. The initial approximation of the pressure 
perturbation required to influence the flow downstream of the 
LEBU was found close enough to the optimum value for enhanced 
large eddy cancellation (it deviated by 2 dB). The optimum 
values of the experimental parameters were very close to 
those initially projected with the concept of acoustic 
excitation of a boundary layer plate manipulator. Thus, for 
further verification, boundary layer velocity profiles were 
measured downstream of the LEBU when the processor operating 
parameters were set at their optimum values. 

3.4. MEASUREMENTS WITH AN OPTIMIZED CONFIGURATION 

The operating parameters of the processor and power 
amplifier were set to the previously established optimum. 
Detailed boundary layer velocity profiles were obtained 
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downstream of the plats manipulator along the midspan of tha 
test surface. Tha momentum thickness was calculated and 
plotted in figure 29 as a function of downstream axial 
distance x. The three configurations compared were 

i. The plain flow 

ii. The acoustically excited LEBU, and 

iii. The optimally excited LEBU. 

The axial variation of momentum thickness exhibited a 
trend similar to the initial configurations. The momentum 

thickness close _to_ the_- LEBU trailing edge for the 

optimally excited flow was again slightly higher than tha 
plain flow's but lower than the simply excited configuration. 
The slope of _the- curve is also smaller- and the resulting 
momentum thicknesses downstream are less with optimum 
acoustic excitation. 

The values of Q were least square curve fitted to a 
power law of the form 0= ax b and the skin friction 
coefficient was estimated using c f =2abx b \ Numerical values 
of the coefficients a and b for various flow configurations 
are shown in table 2. According to these values of the skin 
friction coefficient, optimized acoustic excitation can 
enhance the LEBU's effect by reducing Cf between 2 and 0 per- 
cent with respect to the initially excited configuration and 
by 10 to 17.3 percent when compared to the plain flow 
conf iguration . The improvement of about 2 percent 
additional c f reduction over the unoptimized case was 
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moderate. -This is because the operating parameters of the 
initial configuration were set close to optimum and because 
the least square curve fitting procedure is itself an 
averaging process. 

There may be some experimental error in the measurement 
of Q when the absolute numbers of figures 29, 14 and 15 are 
examined . It - must — be noticed- that always during the data 
acquisition every point of the velocity profile was measured 
with and without acoustics. All the data indicate that 
excitation reduces „ drag and optimized ^excitation improves 
.slightly- _the. effectiveness -of. .the -pulse .-- Some of the data 
could be combined to show results superior to those given in 
figure 29, but the data shown are typical. 
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Table 2a. Curve fitting parameters for the acoustically 
optimized LEBU #5 configuration. 


Configuration 

““plain 

LEBU #5 
( excited) 

LEBU #5 

(optimum acoustic 
excitation ) 


a b 

0.002442 ---0.35175 
0.002912 0.717327 

0.002776 0.731555 


r 

0.993 

0.992 

0.995 


x[ra] 

“ 2.0 

2.5 
- 370 - 

3.5 
4.0 


Table 2b. Axial variation of the skin friction 
coefficient, c^ , for the optimized 


C f (PL) 


0.003754 

0.003532 

0.003535 

0.003455 

0.003337 


configuration. 

,c { (L5A) _ . ... c { ( P L ) - c 4 ( L 5 A ) _ c f (L5 A^ 
Cf(Pb.) 


0.003425 

0.003213 

0.003049 

0.002917 

0.002308 


0 . 037 
0.115 
0.137 
0.155 
0.171 


0.003373 

0.003177 

0.003025 

0.002902 

0.002300 


c,(PL) -c < (L5 App) 

CfCPL) 

0.101 

0.125 

0.140 

0.160 

0.173 


PL signifies the plain flow conf iguration 


Table 2c. Axial variation of the total drag coefficient, 
c D . , for the optimized configuration. 


x[m] 

c^PL) 

c p ( L5A) 

C,(PL)-Cr( L5A) 
cjPL) " 

CpC L 5 A 0 ^) 

Cn(PL)-Cr( L5Anp) 
C^PL) 

2.0 

0.004407 

0.004738 

-0.035 

0.004610 

-0.046 

2.5 

0.004254 

0.004495 

-0.054 

0.004342 

-0.013 

3.0 

0.004150 

0.004259 

-0.029 

0.004134 

0.004 

3.5 

0.004055 

0.004037 

-0.008 

0.003957 

0.022 

4.0 

0.003977 

0.003936 

0.010 

0.003327 

0.037 
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3.5. NET DRAG REDUCTION 

The momentum thickness, Q, at an axial distance x from 
'the'test surface leading edge can be used to calculate the 
total skin friction drag coefficient, c D , at this point [2]. A 
momentum balance is applied considering a control volume 
bounded by the test surface and the undisturbed flow, and 
extending downstream from the test surface leading edge to an 
axial location x. Then the total drag coefficient is given 
by: 

c D (x)=2.0-§!iSl - ' , (3.1) _ 

D x 

”■ - Ti'fgu're' ~29 is^Teconsidered ^here ^observing that the 
momentum thickness corresponding to an acoustically excited 
flow is lower than for the other configurations. According to 
equation 3.1 and the data listed in table 2, a 4 per cent 
net drag reduction, with respect to the plain flow, can be 
obtained at a distance of 50$ downstream of the LEBU 
manipulator. The acoustic power input to the devices was not 
taken into account. 

The actual results of optimization are, however, more 
beneficial than those revealed by table 2. Comparing figure 
15, which is the unoptimized L5 case with figure 29, which is 
the optimized L5 case, two important differences can be 
noted. First, the momentum thickness near the trailing edge 
of the blade is larger with excitation in both cases, but the 
optimized case shows Q only slightly higher with excitation. 


xiooa 
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The more effective large eddy cancellation imparts less 
initial drag into the system when properly optimized. 
Secondly, for the optimized case, net drag reduction occurs 
-earlier- - than -the unoptimized-case and-over a significantly 
larger area. In the non optimized case of figure 15, net drag 
reduction is observed starting at approximately x=3.5 m. For 
the same flow and same configuration but with carefully 
optimized acoustic input, the region where net drag reduction 
begins has moved upstream to approximately x=2.9 m (figure 
29). Measurements- in- the downstream direction beyond the 
• range -o f j 60.-(£ wsreGl im i ted- by-; the -.a vai labie length of the wind 
-.tunnel test section. 

3.6. SKIN FRICTION MEASUREMENTS __0F VARIOUS LEBU 

CONFIGURATIONS 

Prior to the experimental set-up L5 (LEBU #5) various 
large eddy break-up devices were studied in a more 
preliminary manner to establish the effectiveness of a single 
plate and, the effect of an acoustic pulse upon the flow 
downstream of a large eddy manipulator. It is appropriate 
that some of this preliminary data should be presented here. 

3.6.1. PRELIMINARY LARGE EDDY BREAK-UP DEVICES 

All the boundary layer manipulators were mounted at a 
distance of 2.0m from the test section leading edge. The wind 
tunnel velocity was set for a unit Reynolds number of 
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700,000 m 1 and the pressure gradient throughout the flow 
was zero. The height of the boundary layer at the LE3U 
location was approximately 4 cm above the test surface. No 
-.upstream hot-: -film . sensor was installed - and no acoustic 
excitation was applied. The geometry of these manipulators is 
given in table 3. 

Table 3* Geometry of preliminary single LEBU’s 


conf iguration 

height above 
floor [<53 

chord 
length [(53 

thickness 

[$] 

LEBU#1 

- -JT00 

. 25 

0.0038 

LEBU #2 

- '0183 

' “1.215 

0.0033 

LEBU#3 

0.50 

0.95 

0.0050 

Boundary layer 

velocity 

profiles were 

measured 


downstream of the various LEBU plates along the centerline of 
the test surface. The measurements extended a distance of 63 
boundary layer thicknesses from the LEBU, with a spacing for 
each axial station of SS'. The values of the momentum 
thickness, Q, for the plain flow case and the various plate 
manipulators are shown in figure 30. 

The values of Q for all the axial stations were least 
square curve fitted in a power law of the form 0=ax b and 
the local skin friction coefficient was estimated following 
the method of section 2.3. Table 4 contains the numerical 
values of the curve fitting parameters a,b 3nd r 2 for the 
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Figure 30. Momentum thickness, u, versus axial distance 
x, for various preliminary configurations. 
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skin friction drag reduction over a distance of Szfr from 
the LEBU. 


Table 4. 


Curve fitting parameters for 
figurations ( was -measured 
625" from LEBU). 


various flow con- 
to a distance of 


conf iguration a 

PLAIN FLOW 0.002565 


b 

0.715993 


percent drag 
reduction 
at 6 2(5^ 


0.994 


LEBU #1 


0.003005 0.6595 


0.991 -3.0 


. LEBU" #2 0.003104 - - 0.62774 

ITEBU' ?3 003094!-.“-: 0. 62952. r 


0.953 

-0.994 


2.7 

273 


According to the data listed the achieved drag reduction 
was negative or moderate. These measurements were preliminary 
and mostly used to develop the experimental set-up and the 
data acquisition procedure. The configuration LEBU #3 was 
considered as reasonably effective device and it was decided 
to proceed with the acoustic excitation after achieving some 
moderate success with boundary layer manipulation. 

3.6.2. PRELIMINARY ACOUSTIC EXCITATION 


The acoustically excited thin plate configuration was 
designated as LEBU #4. The large eddy detector probe was 
installed upstream of the LEBU at a distance of 2. 4(5'' from 
its trailing edge and at a height of 0. 9^ above the wall. The 
time delay of the processor was set at 9.5 msec and the 



pressure level of the acoustic wave at the LEBU's height was 
100 dB. The geometry of LE3U #4 (chord, height, thickness) were 
identical to LEBU #3 as was the location of device and the 
wind tunnel operating unit Reynolds number. 

Boundary layer velocity profiles for the following 
configurations were measured downstream of the LEBU at 
midspan axial locations over- a - downstream distance of 60<5- 

i. Plain flow (hot film sensor was installed) 

ii. LEBU #4 (hot film sensor installed but not 
- . - • . : operating) 

,-J "" ’.LEBU ^4 acoustically excited. - *' - 

3 ' ~ = ” ~ 1 T he c a 1 a u 1 a t e d - v arTues: x>£~ m crnren turn - thick ness are plotted on 

figure 31. No major changes between the configurations can be 
observed in the rate of momentum thickness growth and the 
local skin friction drag reduction was moderate. The values 
of Q were fitted to the power law and the local skin 
friction coefficient was calculated according to the method 
of section 2.3. The curve fitting parameters and the achieved 
drag reduction are shown on table 5. 


Momentum thickness,0x 1 0 (m) 



Axial distance.x (m) 


Figure 31. Momentim thickness ,Q , versus axial distanc 
x, for preliminary acoustic excitation. 
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Table 5. Curve fitting parameter and per cent drag 
reduction for LEBU #4 configuration. 


Configuration 

a 

b 

r 

percent 

plain flow 

0.00246 

0.94356 

0.992 

drag reduction 

LEBU #4 . 

0.002676 

0.91353 

0.993 

-1.7 to 0.5 

LEBU #4 

0.002699 

0.904223 

0.995 

0.0 to 2.2 

(acoust. excited) 

No net drag reduction was ac 

hieved and 

it was 

decided, based 


rrupon f low. Visualization, _to eaiae the-LEBU to 0.9. - in order 
-“ter b~e't t a r in te ncelp_t = thaL _i.nc.Lde rLt._La r ge eddies. This 
. configuration "was LE3U #5 'which "is 'extensively described in 
the previous sections. 

These preliminary data with little success are the 
result of much time and effort to obtain drag reduction with 
acoustic excitation. This method is not well understood and 
requires a considerable effort to find a configuration which 
works well. Because of inconsistencies with these early 
efforts, it was decided that another independent 
con f i gur a t i o n must be tested to provide additional 
verification of the phenomenon. 
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3.7. ACOUSTIC EXCITATION AT UNIT REYNOLDS NUMBER OF 
10 6 m -1 

To further verify the generality of the large eddy 

cancellation process the wind tunnel velocity was set at a 

6 — 1 

unit Reynolds number of 10 m and LEBU #5 was replaced with 
a different manipulator. LEBU #5 was installed and tested in 
order to examine ^the effectiveness of the acoustic excitation 
at a higher speed. The detailed description of conf iguration 
LEBU #5 is given in table 6 and compared with the 

,i* 

optimized LEBU #5. 


successful 
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Table 6: Experimental parameters for optimized LE3U #5 
and LEBU #6 configurations. 


EXPERIMENTAL PARAMETER 
Unit Reynolds number 


LEBU #5 
6 -1 

0.7x10 m 


LEBU #5 

6 -i 

10 m 


Typical mean flow velocity ... 
Boundary layer sandpaper trip length 
Boundary layer trip location ~ i '* 

LEBU trailing edge location 
LEBU material 

LEBUchord length , 

-LEBU .height -above test surface" 

LEBU thickness - - - - ----- ^ 

Boundary layer thickness at 
LEBU location 

Acoustic input hole diameter 
(hole location at LEBU’s trailing 
edge , midspan) 

Upstream eddy detector distance from 
LEBU’s trailing edge 

Tima delay between large eddy 
detection and acoustic pulse 

Upstream eddy detector height 

Momentum thickness Reynolds number 
at LEBU location 

Expected large eddy 
passing frequency, Hz 

Predicted pressure perturbation at 
LEBU's near wake 


11 m/ s 

15.5 m/ s 

21 cm 

31 cm 

3 to 24 cm 

3 to 34 cm 

2 . m 

2 . m 

steel 

steel 

<—0/93 S' 

1.04(5' 

:~77q£bo& 

0 .75$.. 

z -.0:0050 S' 

0 . 0054 s 

~4.1 cm 

~3.7 cm 

3/3 inches 

3/9 inches 

2.A5S 

2.5 S 

10.5 msec 

5.5 msec 

0.85^ 

o.sScJT 

3,100 

3,700 

110 

1 65 

100 dB 

104 d3 


SPL of acoustic pulse at LEBU’s 
height (2 dB higher than predicted 
pressure perturbation). 


102 dB 


105 d3 
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Turbulent boundary layer velocity profiles were measured 
at midspan locations downstream of the LEBU manipulator as 
usual. The momentum thickness, 0, was calculated as described 
in section 2.3 and plotted versus downstream distance, x, from 
the test surface leading edge. Figure 32 shows the axial 
variation of Q for the usual three configurations. 

The momentum thickness of the manipulated flow is 
greater than the plain flow case because of momentum loss 
' G * r '~ T ~ imposed bjr the- embedded thin plate. -The ^-acoustically excited 

'Case -produces, u-jowar. m-araTentuml loss! at- this location. This 
: v-vc indicates “that the : acoustic “pulse does - not subtract kinetic 
energy from the flow and it improves the efficiency of the 
LEBU conf iguration in the wake region. From figure 32 it can 
be observed that the momentum thickness of the acoustically 
excited conf iguration again grows at a slower rate compared 
to the simply manipulated boundary layer. Further downstream 
at x = 3.75 m (47^ from the LEBU) the momentum thickness of 
the acoustically excited configuration falls below that for 
the other conf igurations indicating net drag reduction. 

The variation of the skin friction coefficient, c^ , with 
the axial distance from the test surface leading edge is 
shown in figure 33. All measurements were obtained at the 
test surface raid span aligned with the upstream sensor and 
the acoustic wave port. The numerical values of the skin 
friction coefficient were calculated according to the 




Figure 32. Momentum thickness, □, versus axial distance 
x, for various flow configurations at unit 




nx . c = 2abx 
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procedure described in section 2.3. All least square curve 
fits again had a correlation coefficient above 0.99. The LEBU 
configuration and the acoustically excited flow were compared 

. - to _ the., plain _ flow ( hot film sensor installed ). According to 

the numerical values of the skin friction coefficient, 
acoustic excitation can enhance the LEBU's effect by reducing 
... _..Cf between 3.3 and 6.0 percent with_respect to the LEBU 

configuration. Comparison to the plain flow shows a reduction 
of the wall shear stress between 6.5 and 13 par cent. In 
->x'^£5-?3dditionr..qa 2. percent net drag-reduction= was measured at a. 
dlstanceF."of . 54 $ from the LE3U. f 

. the R 7.3 Aithougha-i: the: data, arse not -..as: impressive as the previous 

’ • A • optimized co.nf iguration , they do suggest that the mechanism 

of enhanced large eddy cancellation is real. Anders [13], has 
shown that variations of the manipulator microgeometry 
exercises a strong influence upon the results so absolute 
numbers are not extremely meaningful. The trends of the data 
are, however, consistent with previous findings. 

Space-time cross correlation functions (CCF) were 
obtained using the hot film eddy detector probe upstream of 
the LEBU plate as a fixed reference location and another 
probe at a distance of 4.0,^ boundary layer thicknesses 
downstream (figure 17). The height of the downstream sensor 
was 0.9^. Figure 34 shows representative space-time cross 
correlation functions for the plate manipulated and the 
acoustically excited boundary layers. When 


acoustic 


0.084 
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excitation is applied the downstream flow is again correlated 

to the upstream point signif icantly less. The peaks of the 

CCF correspond to eddy convection time from the upstream to 

the downstream sensor. This indicates that acoustic 

excitation enhances the destruction of the large eddies which 

are sensed at the upstream location. These measurements are 

in good agreement with similar results of configuration LE3U 

tf 5. In addition measurements of the low frequency turbulence 

spectrum, downstream of the plate manipulator showed that 

iacousticrrexcitation can=ceduce the RMS^value of the axial 

.velocity fluctuation approximately 2 percent. 

: : z Figure. ~35 “shows ~the -variation of— the ratio r _ = ^ 

- - ° 0(0 

with the spanwise distance from the test surface 
ce n te r 1 ine , z , at a distance of x=1.76m from the LEBU. 
According to previous measurements of momentum thickness, Q, 
(figure 31) the ratio r^ was found smaller than unity at this 
axial station. In this location it was found that the 
acoustic excitation persisted with a tendency to spread 
spanwise at a half angle of about 1.2 degrees. The reduction 
of the angle of acoustic influence when compared to LEBU #5 
is attributed to the increased mixing process of the 
turbulent boundary layer due to higher free stream velocities 
in which LEBU §S was tested. 

Proper acoustic excitation of a thin flat plate 
operating in a higher unit Reynolds number flow was again 
found effective in enhancing the large eddy break-up process. 



0(with acoustics ) 
“9 (without) ’ 
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Figure 35. Spanwise variation of momentum thickness ,Q, 
for the excited and unexcited LEBU at unit 
Re = 10 6 m" T . 
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The processor operating parameters were set according to 
higher free stream velocities. The achieved skin friction 
drag reduction, reduction of flow correlation across the 
LEBU, reduction of the low frequency turbulence and the 
moderate net drag reduction were consistent _with the results 
obtained from LEBU#5. This reinforces the initial concept 
that proper acoustic excitation of a plate boundary layer 
manipulator can lead to enhanced large eddy destruction. 
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4. CONCLUSIONS AND RECOMMENDATIONS 

4.1. CONCLUDING REMARKS 

Large eddy structures in-a turbulent boundary layer can 
be effectively manipulated by acoustic waves resulting in a 
net reduction in skin friction drag downstream of a 
.manipulating plate.. The. vortex unwinding, at the LEEU's 
trailing edge when the blade encounters a large eddy is 
assisted by acoustic waves. These acoustic pulses must be 
.phase unlocked ^toitheieonveetedteddies ^arrive at the LEBU - 
"trailing >s edge ^together)' and -must-have '-.an . adequate sound - 
pressure level- in border -to' be effective. The coherent 
■structure : _was .rassumed^to be ‘related itoi the Reynolds stress 
production and the sublayer bursting “events. Thus, The 
additional elimination of the large scale structure caused by 
the acoustic waves improved the effectiveness of a single 
plate manipulator. Skin friction drag reduction was achieved 
for two different low Reynolds number LEBU configurations. 

The skin friction coefficient c f was reduced as much as 18 
percent at a distance of 50^ from the LEBU when the 
acoustically excited plate configuration is compared to the 
plain flow. In addition, a moderate net drag reduction was 
achieved at about the same distance from the manipulator. 

The application of the proper acoustic input made the 
flow less correlated across the LEBU, and reduced the level 
of the low frequency part of the turbulence spectrum of the 
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manipulated boundary layer. In addition, flow visualization 
of the acoustically excited boundary layer showed the further 
reduction in coherent motion. The initial hypothesis, that a 
phase _l.ac-k.ed- -.acoustic input to the flow can reduce the 
turbulence mixing and the wall skin friction was verified. 

The effects of the acoustic excitation spread slowly in 
-the - span.wisa_direct ion as the flow: jconve-Oted downstream. This 
was verified by carefully measuring the momentum thickness, Q, 
in the spanwise direction downstream of the manipulator. The 

i ™ * r„is pra a d Lng wa s veri f ted^wi t h r^eiPordity: carccs s.cor. relation -a n d 

turbulence ' "-intensity: measure are rrt s'. v^T he "acoustic excitation - 
i -effect ‘spreads .at Va -half angle -of 2 degrees for the most 
--- effective configuration.; The "Spreading angle was reduced as 

- - Re -and turbulent mixing increased. It was assumed that 
secondary flow vortices induced by the acoustic pulse are 
convected into the turbulent boundary layer and cancel some 
of the "wsaksned" large eddies passing the plate manipulator. 

An alternative explanation for the spreading process can be 
based upon the assumption that the large scale structure and 
the sublayer bursting phenomena are interrelated and preserve 
each other. The passing of a large eddy causes sideways 
eruptions of fluid. These bursting and sweeping motions cause 
the generation of the large scale structure. The interruption 
of this chain of events by eliminating the large eddies at a 
spanwise position would yield has a spreading effect upon the 
reduction of bursting events and further generation of large 
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The operating parameters of the processor device which 
controls the acoustic input could be set to an optimum in 
order_to_scc_ount for the— large-. eddy.- convection velocity,, 
their particular frequency of appearance and the level of 
the velocity fluctuation they represent. The optimization of 
the processor operating parameter was based upon the observed 
phenomenon of reduced correlation across the LEBU when 
acoustics is applied. For maximum large eddy cancellation 
=tha ; ur_eferencs- yoltage .mb us£ r.psa.t; atftipr oduce. approximately 

the~,S3tne~ numb-e-r- o-f- -accru stic' -pulses ; per unit time -as the 
number xsf predicted large, eddies pec cuTidt time. The process 
of.-.reduced correlation is very much dependent upon the 
setting of the reference voltage and must be set within 10 
millivolts from the optimum value for best results. The flow 
is least correlated when the time delay is set equal to the 
large eddy convection time from the detector probe to the 
LEBU trailing edge. Also, the initial approximation of the 
pressure perturbation required to influence the flow 
downstream of the LEBU was found close enough to the optimum 
value to achieve enhanced large eddy cancellation. The 
optimum values of the experimental parameters were very close 
to those initially projected. These results verify the 
concept of increased drag reduction by acoustic excitation of 
a boundary layer plate manipulator. 


4.2. RECOMMENDATIONS FOR FURTHER INVESTIGATION 


Acoustic excitation was proven effective in improving 
the plate manipulator effect upon a turbulent boundary layer. 
The basic concept of the large eddy cancellation by acoustic 
waves should be studied in the case of a manipulator plate in 
inviscid flow. A LEBU device should be positioned in 
undisturbed flow so that it may encounter vortices 

artificially generated upstream. Acoustic pulses must be made 
to arrive at the blade trailing edge together With these 
-vortices ]• ^---The . -effect - of -acoustic waves —upon the 
iresulting . flow. field .could .be examined^to further establish 
the basic mechanism -developed in this study. 

Further measurements, ...... beyond the_range of 50^ are 

recommended in order to establish the extent of the 

persistence of the acoustic effect downstream of the LEBU. 
Extensive measurements of the axial and longitundinal 

velocity fluctuating components u' and v* should be taken to 
determine changes in the turbulence energy production. The 
production of skin friction drag is closely related to the 
Reynolds stress component -pu'v'; this quantity should be 
also examined comparing an acoustically excited LEBU 
configuration to single and tandem plate manipulated 
turbulent boundary layers. 

Further investigation of the spreading tendency of the 
acoustic effect upon the flow regime downstream of the LEBU 
should be pursued. Two hot film "scout" probes can be 
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installed upstream of the plate and the acoustic pulses can 
be input from two different ports, located relatively close 
to each other. Boundary layer traverses will show if both 
zones of acoustic influence merger downstream of the LEBU. 

A LEBU with two plates in tandem^exposed to acoustic 
- waves - may be also examined because this configuration has 
been found most promising without excitation. The boundary 
layer upon a LEBU plate is laminar. Possible separation of 

.• — the _.LE£U. Jaouodary.. layer., can. be. -.prevented by forcing it 

- , . * ; -..i- . j - ± u rbii le n t. ,wLth,a -tr,iph wire .at.. the thin* -plate leading edge "or ' : 
-■r-i ' :g i - hy Toughening its" surface . - Such a LEBU configuration is also 

idU .worth studying when exposed to acoustic pulses. 

— — r -z Most' cfr.'--tha: previous experiments with plate manipulators 

were carried out at a slightly higher speed than the present 
investigation. An experimental setting at higher velocities 
will require an acoustic pulse of higher frequency and sound 
pressure level. Improvements and redesign of the pressure 
pulse mechanism must be undertaken. Such a configuration, 
closer to real flow conditions, in conjunction with a 
variable pressure gradient and airfoil type LEBUs should also 


be studied. 
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6. APPENDICES 


6.1. FLUCTUATING PRESSURE AT THE LEBU’S TRAILING EDGE 


Based upon the procedure suggested in reference [32], a 
prediction of the fluctuating pressure at the LEBU’s trailing 
edge was obtained. ~ This- permitted the time averaged sound 
pressure, level from 0 to 2000 Hertz of the acoustic waves 
impinging upon the flat plate to be set at a value which 
could influence the flow. The study of reference 32 examines 
a flat plate in an inviscid flow. The resulting equations are 
based “"-upon-- the solution .::of--- the --complete Navier-Stokes 
Equations using the Theory , of .Asymptotic Expansions., . 

:1 . The nondimensional perturbation pressure gradient at the 
LEBU trailing edge_ was found as (fig. 7, reference 32): 

P(TE)=-0.301 (5.1) 


The above quantity is defined by: 
P(T£)= —a D ? 2 . 


(5.2) 

y z'Uo 

where, DP is the fluctuating^ pressure . The parameter e for 
laminar flow is equal to Re c 8 , where, Re c is the Reynolds 
number based upon the LEBU's chord length. The constant 
C =0.33205 appears in the Blassius solution for a flat plate. 
The flow density, p , was given a typical value of 1.25 Kg/m 
and U 0 is the free stream velocity. Considering LEBU #5 
mounted at a height hsO.ScJT, U 0 was approximated as 10.5 m/s, 
and the parameter -^2. was 670,000 m at the LEBU's height. 

The Reynolds number, Re c was found to be ~25,000 and 
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the boundary layer upon the LEBU was assumed to be laminar. 
The transition Reynolds number Re x , with a 3 per cent 
turbulence intensity (typical at LEBU's height ) was 
approximately 75,000. 

After substituting the given quantities into (6.2), DP 

■ - was z- found equal r -to 99.95 decibels. -Similar calculations for 

- . . the LEBU #5 resulted in a pressure perturbation of 104.3 dB. 

A sound wave with such a low frequency ( 0 to 200 Hertz ) 
time averaged pressure level was considered adequate to 

— -.^trigger . the shedding of a vortsx-at the_LEBU's trailing edge 

„ _ _ ~~a ri d-i n f 1 ue h c el a rge^ed d y cancellation. 
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6.2. PREDICTION OF TURBULENT BOUNDARY LAYER. 
6.2.1. INPUT DATA FOR HEAD'S METHOD 


Head’s method [23] was utilized to predict the growth of 
characteristic lengths of a plain turbulent boundary layer. 
Head’s method requires the following data as input: 

_1, The momentum thickne_ss ~a t t he beginning of the 

predicted flow regime (end of sandpaper trip). 

2. The initial shape parameter, H. 

Aloe Nfr'ee streams velocity, U 0 as a function -of 
LClm ‘.The flow kinematic viscosity assumed constant for 


this case. 

According to “equation- (3.1 )- the total drag coefficient c D 
at the end of the sandpaper trip is defined as: 


C n (x)= 2 


0(X) 


(5.3) 


’D v "' - X 

The drag coefficient at the end of the sandpaper trip was 


empirically predicted [1] by: 


• Re Y 


0.4 A. 


0.3 Re x ^ x ) 

A 




A 


) 


(6.4) 


In the above relation Re^ is the Reynolds number based 
upon the sandpaper trip length x, K is the average roughness 
height of the sandpaper and \ is a parameter defined as 


\-c 

The skin 
sandpaper 


-£ — )°' 5 

c fW 

friction coefficient c f at the 
trip was predicted by the empirical 


end of the 
formula from 


White [ 1 ] . 

Y -02 

c f (x) = (l .4+3.7--) 


(5.5) 


SpSffAI 

0!? POOR; 


* page is 
QUALITY] 


The average sand-paper roughness height K was about 0.001 
m, the sandpaper trip length was 0 . 24 m, and the wind tunnel 
flow velocity unit Reynolds number _ was 700,000 m"' . Proper 
manipulation of the above "formulas can provide a prediction 
of the momentum thickness at the end of the sandpaper trip. 

~ The initial or starting .values of H in Head's method are 

re lat:ivai.y.-unri m po r t anir-rex cep t — f o r t he first few points of the 

calculation. A shape parameter for a flat plate flow H=1.375 
was assumed by averaging measured values of H along the flat 
— J '"plate-.~'The -free -stream -Vertoc irty “was ^considered to be 11 r m/ sr 
'■-'•The growth of the momentum th ickne s s , 0 ipredicted "By Head’ s 

' method . fo.c ia ...plaia__ flo.w configuration is compared with 

- - measured s in figure 9 . 


6.2.2. VELOCITY PROFILE OF A TURBULENT BOUNDARY LAYER 

The logarithmic and the outer region of a plain 
turbulent boundary layer were predicted using a modified 
Coles velocity profile [33]. The velocity u at a boundary 
layer height y is thus given by: 

^ = X Un(J 7 !l, + ( T' )2(1 --7- ) + rrw ' ! /- >1 * B ’ (5 ' 6) 

where, u t is the wall skin friction velocity, IT is the wake 

y 

parameter and w,( ^— ) is the well known wake function. The 
logarithmic velocity law constants k and B were assumed equal 
to 0.41 and 5.0 respectively. The wake function [ 34 ] is 

defined as: 



( 6 . 7 ) 
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which also can be approximated with the sine function 
proposed by Coles [1]: 

w(^-)=2sin 2 (yy) . (6. S) 

The raore^ complicated relation (6.7) was used because it 
g3ve better agreement with the measurements. 

At the wall the relation (6.7) yields w(0)=0.0,and at 
-the .'redge--of the-boundary- layer . .0)=2.0, Relation 

(5.6) then yields: 

~ n-2.).B (6.9) 

U t k V 

Subtracting-equatlon--(6?6-) from (6.9) one can get: ” r 

t Uolilrw ^tlri(^-(^j 2 (V*-^)*TT(2.-w,(^-))] (6.19) * - 

...A--. The., wake- parameter. . for.-an-equilibrium boundary layer 
remains constant ; for a flat plate measurements suggest TT= 0.55 
[1]. The wall skin friction velocity u t was determined from 
the definition rearranged as: 

u t =Uo ("£* ) (5.11 ) 


The local skin friction coefficient 
formula [1 ] 


c, = 


0.3 e 


-1.33H 


'■(logRe/ 7 ** 0 ' 3 ^ 

The boundary layer velocity profile 


was predicted by Coles 

( 6 . 12 ) 

predicted by using these 


relations in equation (6.10) is compared with a measured 
plain flow profile in figure 11. 
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' 6 . 2 . 3 . BOUNDARY LAYER VELOCITY PROFILE FOR A 
CONFIGURATION 


LEBU 


u LE =u pL - du w , 


« d .V5. du ma« 




The velocity, u L ^, at a height y of the LEBU 
manipulated turbulent boundary layer was predicted as: 

(6.13) 

-where, _u pL is the velocity predicted for a plain configuration 
and du w is the velocity defect caused by the LEBU plate wake. 
For inviscid free stream flow and a two dimensional wake, du w 
is .estimated -by the relation [2]: 

~ (6.14) 

where , du ma ^'is;-the ^maximum “velocity idefec-t at the wake's, axis 

of ^symmetry ,: ^3( X ) is the wake serai-width and y A ‘ is the 

zlongitundinal - distance from the -wake ’s axis of symmetry. 

According to previous studies [35] for a similar two- 

dimensional wake the quantities ft and dUj^^at a distance x 

from the plate trailing edge, are defined as follows: 

A*) . x 0.5 

CTT =( ^ ) (5.15) 

6(o) ©(o) 

and 

du — v -0.5 

) , 


= (- 


( 6 . 16 ) 


U 0 ’ 0(o) 

where U 0 is the boundary layer velocity at the LEBU height 
(free steam velocity) and 0(0) is the momentum thickness at 
LEBU trailing edge. The quantity 0(0) was considered as 
being twice the momentum thickness of the laminar boundary 
layer grown upon the LEBU.' The Blassius solution [2] for a 


flat plate with a laminar boundary layer suggests 

0(c)-O.664 )° ’! 

u 0 


(5.17) 
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In the above relation c represents the LE3U chord length, and 
the parameter — ^2- was 670,000 m" 1 (at the LEBU's height). 
-The logarithmic velocity law constant k was assumed to retain 
the value of 0.41 [36], and the wall skin friction velocity 
was calculated frpra__(5, 10) . The skin -friction coefficient 
resulted from previous measurements (section 2.3). A 
"theoretically - predicted - velocity profile for LEBU #2 is 
compared with a measured profile on figure 13. 


:r 4v PETER MIRATION, QF--e- f nFROM -VELOCITY- PROFILE 


-mr--'The-’ following me thod. .of calculating; ;Cj is a variation of 
the procedure -suggested "by Bradshaw £24];.. In the logarithmic 
-region of a boundary layer, the. velocity u at a distance y 
from the wall is given by: 

(5.13) 

where u t is the wall friction velocity , v is the kinematic 
viscosity and k and B are constants equal to 0.41 and 5.0 
respectively. For a flat plate the logarithmic region 
extends approximately over 33<-^r 1 <350. So, for 


u i • r y u t^ 0 
= 2L_ ln(— t) +B ; 


-^ = 200 

relation (6.18) is valid and yields 
-^-=17.923. 

u t 

The skin friction coefficient is given by: 

— -(-^) 2 
2 1 U 0 * 

where U Q is the free stream velocity. Only one point (y,u) on 


(6.19) 


( 6 . 20 ) 


( 6 . 21 ) 
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the logarithmic region of a measured boundary layer satisfies 
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all the above relations. Determining this point can provide a 
value for at a certain axial location. After traversing 
across a boundary layer various values of u are substituted 
intxr-re-l-a-ticm— (S.SO^r'The-x^esul'ting u and corresponding y 
are used to determine the point on the logarithmic region 
which gives closest to 200. This requires a local least 

-square fitting of the form 

u=ay b (5.22) 

which is applied for measured values of u and y in the 
vicinity of _this point.- T h us ; though 'the pr ac i s'e ' po in t 
where y*.=200.. mayilriot have teenLjiieLasured.u the corresponding u 
and y for ^^ = 200 can -be_-approx imated very accurately. 


Manipulation of ithase formulas yields: 

y=( MMLL.)Ub (5.23) 

o 

The resulting value of y and u satisfies the relation (6.20) 
and can be used to determine c^ from (6.21). 

Example calculation: 


Measured quantities at x=2.75 m, z=0.0 m. 
-^=700,805 m“1 u =1.614 10~ 5 , U 0 -11.311 m/s 


y[m] 

u[ m/ s] 

u^m/s] 

yu t 

0.005 

7.529 

0.420 

130.1 

0.005 

7.700 

0.430 

159.703 

0.009 

8.033 

0.443 

222.376 


Least square curve fitting of the measured values of u at a 
distance 0.004 to 0.016 ra from the wall results in 
u=0.007256y^' 1556 with correlation coefficient equal to 
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0.993. Then substituting into the proper relations one 
get: y=0. 007266 m, u=7- 96l056m/s, u t =0. 4442m/s, and 
c f =0.003084. The value of c f is the desired result. 


can 
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